Statistical nature of cluster emission in nuclear liquid phase 
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The emission of nuclear clusters is investigated within the framework of isospin dependent lattice 
gas model and classical molecular dynamics model. It is found that the emission of individual cluster 
which is heavier than proton is almost Poissonian except near the liquid gas phase transition point 
and the thermal scaling is observed by the linear Arrhenius plots which is made from the average 
multiplicity of each cluster versus the inverse of temperature in the nuclear liquid phase. It indicates 
of a statistical nature of such cluster emission in the models. The "emission barriers" which are the 
slopes of the Arrhenius plots are extracted as a function of the mass or charge number and fitted 
by the formula embodied with the contributions of the surface energy and Coulomb interaction. 
The possible influences of the source size, Coulomb interaction and "freeze-out" density and related 
physical implications are discussed. 

PACS numbers: 25.70.Pq, 05.70.Jk, 24.10.Pa, 02.70.Ns 
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I. INTRODUCTION 

Hot excited nuclei with moderate temperature can 
be formed by the collisions between heavy ions at low- 
intermediate energy and they finally de-excite by the dif- 
ferent decay modes, such as the light particle evaporation 
and the emission of multiple intermediate mass fragment 
(IMF), ie. multifragmentation. Even though extensive 
studies on multifragmentation have been carried out ex- 
perimentally and theoretically, it is still an open ques- 
tion to clarify definitely whether the multifragmentation 
is statistical or dynamical, sequential or simultaneous. 
Among such efforts, Moretto et al. found that there ex- 
ists the resilient reducibility and thermal scaling in mul- 
tiple fragment emission process, which seems to show one 
a possible interpretation picture to look and understand 
the multifragmentation. Originally, they observed that 
the experimental Z-integrated fragment multiplicity dis- 
tributions P™ are binomially distributed, ie. 
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in each transverse energy (E t ) window, where n is the 
number of emitted fragments and m is interpreted as the 
number of times the system tries to emit a fragment. The 
probability of emitting n fragments can be reduced to a 
single-particle emission probability p which gives linear 
Arrhenius plots (ie. excitation functions) when ln(l/p) 
is plotted vs l/y/(E t ). By assuming a linear relation- 
ship between y/(E t ) and temperature T, they proposed 
that the linearity of the observed ln(l/p) vs l/y / (E t ) 
can be interpreted as a thermal scaling of the multi- 
fragment process |, |, |]. In this case, these linear 
Arrhenius plots suggest that p has the Boltzman form 
p oc exp~ B / T with a common fragment barrier B. How- 
ever, since the binomial decomposition has been per- 
formed on the ^-integrated multiplicities, typically as- 



sociated with 3 < Z < 20, the Arrhenius plot generated 
with the resulting one fragment probability p is an aver- 
age over a range of Z values. Some comments and criti- 
cisms on the above binomial distribution and the thermal 



scaling were also raised 



[s} |(| fjj. One paper, that of 
Tsang and Danielewicz [4], found that the assumption 
of T a y/(Et) may be only valid for compound nuclei 
formed at low-to-moderate temperatures and it will be 
broke down for the intermediate energy heavy ion colli- 

Au at E/A = 37 to 110 MeV 



sions such as 40 Ar 
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due to the mid-rapidity particle emission from the over- 
lap region of projectile and target as well as the delay 
emission from projectile-like and target-like residues. In 
this case the above interpretation of thermal scaling is 
not vaild in the point of experimental view. Moreover, 
they found that the fit with the binomial distribution 
can be replaced by the Poisson distribution in the con- 
straint of charge conservation. As a consequence, the 
reducibility of fragment emission to bionomial distribu- 
tions does not imply the fundamental significance for the 
parameters of m and p thereby extracted. If p is not an 
elementary emission probability, then the Arrhenius law 
will fail also. 

Later, instead of analyzing for Z-integrated multiplic- 
ities, Beaulieu and Moretto et al. analyzed the behavior 
of individual fragment species of a given Z for higher reso- 
lution experimental data and noticed that the n-fragment 
multiplicities P(n) obey a nearly Poisson distribution, 



P(n) = 



< n >' 
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where n is the number of fragments of a given Z and the 
average value < n > is a function of the total transverse 
energy E t , and were thus reducible to a single- fragment 
probability proportional to the average value < n > for 
each Z j8| pJ|. Similarly the < n > is found to be pro- 
portional to exp~ B l T providing that T oc ^jE t , ie. there 
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exists also a thermal scaling law. More recently, Elliott 
and Moretto et al. discovered that the common features 
of Poissonian reducibility and thermal scaling can also 
be revealed in percolation and the Fisher droplet model 

0,0- 

In the present work, we would like to make a theo- 
retical reexamination on the Poissonian reducibility and 
its thermal scaling rather than the bionomial reducibil- 
ity and its thermal scaling. Unlike in experiment, we 
will adopt the true temperature to check the Poisson re- 
ducibility and thermal scaling in the frameworks of the 
isospin dependent lattice gas model (I-LGM) [jl2|, [ll| 
and followed by the isospin dependent classical molec- 
ular dynamics (I-CMD) of Das Gupta and Pan |l4j. Of 
course, we should keep in mind that the assumption of 
\/{Et) a T in experiments may be only valid for com- 
pound nuclei formed at low-to-moderate temperatures, 
but fail in the experimental data at higher temperatures 
In this context, we checked the relationship between 
y/{E t ) and T in the models and found that the assump- 
tion of \f{E t ) oc T is valid at low-to-moderate temper- 
atures where we will study in the following section. By 
investigating the variances and average multiplicities of 
cluster multiplicity distributions as a function of temper- 
ature, we will illustrate that the Poissonian reducibility 
and its thermal scaling is valid for fragment emission in 
the nuclear liquid phase in the framework of the above 
thermal equilibrium models. It may indicate of a statis- 
tical nature of cluster emission. 

The paper is organized as follows. Firstly, we introduce 
the models of I-LGM and I-CMD in Sec. II; In Sec.III, 
the results and discussions are presented. We first show 
some results to support that there exists a Poisson re- 
ducility in the cluster production away from the liquid 
gas phase transition by investigating the ratio of the dis- 
persion of multiplicity distribution to its mean value and 
the Poisson fit to the multiplicity distribution of indi- 
vidual clusters. Second we plot the Arrhenius-type plots 
and find the thermal scaling is valid in the nuclear liquid 
phase. Further we extract the "emission barrier" sorted 
by the different mass number, or light isotope, or charge 
number in I-LGM and I-CMD and use the formula em- 
bodied with the contributions of the surface energy and 
Coulomb interaction to make a systematic fit. The de- 
pendences of the model, the source size, Coulomb inter- 
action and "freeze-out" density are presented. Finally, a 
conclusion is reached in Sec. IV. 



II. DESCRPITION OF MODELS 

Orginally, the lattice gas model was developed to de- 
scribe the liquid-gas phase transition for atomic system 
by Lee and Yang [jl5| . The same model has already been 
applied to nuclear physics for isospin symmetrical sys- 
tems in the grand canonical ensemble ||l6f w ith a sam- 
pling of the canonical ensemble g|, [l| || ||, 
El|, and also for isospin asymmetrical nuclear matter in 



the mean field approximation p2| . In this work, we will 
adpot the lattice gas model which was developed by Das 
Gupta et al. |l2|, pH) . In addition, a classical molecular 
dynamical model of Das Gupta et al. |TJ] is also used to 
compare with the results of lattice gas model. For com- 
pleteness of the paper, here we make a brief description 
for the models. 

In the lattice gas model, A (= N + Z) nucleons with 
an occupation number Sj which is defined Si — 1 (-1) for 
a proton (neutron) or = for a vacancy, are placed 
on the L sites of lattice. Nucleons in the nearest neigh- 
boring sites have interaction with an energy e SiS .. The 
hamiltonian is written as 



A p2 

^ 2m 



i<j 



(3) 



where Pi is the momentum of the nucleon and m is its 
mass. The interaction constant e SlS is chosen to be 
isospin dependent and be fixed to reproduce the bind- 
ing energy of the nuclei [H : 



O.MeV, 



= -5.33MeV. 
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Three-dimension cubic lattice with L sites is used which 
results in pf — j-po of an assumed freeze-out density of 
disassembling system, in which po is the normal nuclear 
density. The disassembly of the system is to be calcu- 
lated at pf, beyond which nucleons are too far apart 
to interact. Nucleons are put into lattice by Monte 
Carlo Metropolis sampling. Once the nucleons have been 
placed we also ascribe to each of them a momentum by 
Monte Carlo samplings of Maxwcll-Boltzmann distribu- 
tion. 

Once this is done the I-LGM immediately gives the 
cluster distribution using the rule that two nucleons are 
part of the same cluster if 



P?/2p 
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where P r is the relative momentum of two nucleons and 
p is their reduced mass. This prescription is evidenced 
to be similar to the Coniglio-Klein's prescription [2a| in 
condensed matter physics and be valid in I-LGM Jl2[ [14 
ppf . To calculate clusters using I-CMD we propagate 
the particles from the initial configuration for a long time 
under the influence of the chosen force. The form of the 
force is chosen to be also isospin dependent in order to 
compare with the results of I-LGM. The potential for 
unlike nucleons is 



« np (r)(-<o) = A \b(—) p - (— ) 9 1 exp(^ — ), 
r L r r i — - a 
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In the above, r = 1.842/m is the distance between the 
centers of two adjacent cubes. The parameters of the po- 
tentials are p = 2, q = 1, a = 1.3, B — 0.924, and A = 
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1966 MeV. With these parameters the potential is min- 
imum at r with the value -5.33 MeV, is zero when the 
nucleons are more than 1.3ro apart and becomes stronger 
repulsive when r is significantly less than ro . The poten- 
tial for like nucleons is written as 

v pp( r )( r < r o) = v np (r) - v np (r ), 
v pp (r)(r>r ) = 0. (7) 

This means there is a repulsive core which goes to zero 
at ro and is zero afterwards. It is consistent with the 
fact that we do not put two like nucleons in the same 
cube. The system evolves for a long time from the initial 
configuration obtained by the lattice gas model under the 
influence of the above potential. At asymptotic times the 
clusters are easily recognized. The cluster distribution 
and the quantities based on it in the two models can now 
be compared. In the case of proton-proton interactions, 
the Coulomb interaction can also be added separately 
and compared with the cases where the Coulomb effects 
are ignored. 

III. RESULTS AND DISCUSSIONS 
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In this paper we choose the medium size nuclei 129 Xe 
as a main example to analyze the behavior of individ- 
ual fragment emission during nuclear disassembly with 
the helps of I-LGM and I-CMD. In addition, the systems 
with A sys = 80 (Z sys = 33) and 274 (Z sys = 114) are 
also studied to investigate the possible source size depen- 
dence. In most case, pf is chosen to be about 0.38 poi 
since the experimental data can be best fitted by pj be- 
tween 0.3po and 0.4p in the previous LGM calculations 
fL2| p4[ , which corresponds to 7 3 cubic lattice is used for 
Xe, (F for A sys = 80 and 9 3 for A sys = 274 . In the 
condition of the fixed freeze-out density, the only input 
parameter of the models is the temperature T. In the 
I-LGM case, pf can be thought as the freeze-out den- 
sity but in the I-CMD case p/ is, strictly speaking, not a 
freeze-out density but merely defines the starting point 
for time evolution. However since classical evolution of a 
many particle system is entirely deterministic, the initial- 
ization does have in it all the information of the asymp- 
totic cluster distribution, we will continue to call pf as 
the freeze-out density. 1000 events are simulated for each 
T which ensures enough statistics. 

A. Poissonian Reducibility 

One of the basic characters of the Poisson distribution 
Eq.(2) is the ratio c 2 ./ < rt,; >— > 1 where a 2 , is the 
variance of the multiplicity distribution and < n.; > is 
the mean value of the multiplicity distribution. The first 
step we are showing is this ratio. We obtain these ratios 
for clusters classified with different masses (^4) , light iso- 
topes (ISO) and atomic numbers (Z) for the disassembly 



FIG. 1: The ratio of ai 2 / < rii > for the clusters classified with 
mass, light isotope mass and atomic number as a function of tem- 
perature. The left panels are for the I-LGM calculation and the 
right for I-CMD with Coulomb. The symbols are illustrated on the 
figure. 



of 129 Xe as a function of temperature in the framework 
of I-LGM and I-CMD with Coulomb in Fig. @. Obvi- 
ously, most of the ratios are close to one, which indicates 
that basically these cluster production obeys the Poisson 
distributions, ie. a cluster is formed independently from 
one another. Of course, we also notice that the values of 
protons are almost lower than the unique, ie. it is nar- 
rower than the Poisson distribution. This could be due to 
protons can be easily produced by an "evaporation" -like 
mechanism in the models, ie. protons can be easily sepa- 
rated from some unstable multi-nucleon clusters. In this 
case, this kind of proton production is obviously related 
to the parent cluster and then a narrower distribution of 
mixed protons could reveal. On the contrary, some points 
are slightly larger than 1 and this behavior becomes ob- 
vious in the mediate temperature range, which could be 
related to the onset of phase transition in such moderate 
thermal excitation |^l[, where the large critical fluctua- 
tion |2f| makes the Poisson reducibility broke-down . 

Besides the above Poisson condition is basically sus- 
tained, the excitation function of the average multiplic- 
ity of rt-multiple individual cluster emission can be well 
fited with the Poisson distribution. For some examples, 
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FIG. 2: The excitation functions of n-multiple cluster emission 
probability (Pn) for elements with Z > 2 emission from the source 
12 ®Xe in the I-LGM calculation. The lines show the expected val- 
ues in different temperature with the Poisson assumption according 
to Eq.(2). 

Fig. H shows the quality of the Poisson fits to the average 
multiplicity of n-multiple individual cluster emission in 
the different temperature for 129 Xe in the I-LGM case. 
In each panel of this figure, we first plot the probability 
of n-multiple emission species ( P(n) ) as a function of 
temperature (as shown by the open symbols), and then 
we connect the Poisson probabilities in the different tem- 
perature as lines in terms of Eq.(2) due to we know n 
and its average value < n > over all n-multiple emis- 
sion in each temperature (as shown by the lines). Obvi- 
ously, these Poisson fits are quite good for almost Z > 2 
over the entire range of T. The similar good Poisson fit 
is overall obtained in the cases of I-CMD. As a conse- 
quence, we think that Poissonian rcducibility is valid in 
the thermal-equilibrium lattice gas model or molecular 
dynamics, which illustrates that the cluster production 
is almost independent each other in the studied models. 



B. Thermal Scaling 

Naturally, we want to know if there exists a thermal 
scaling law in the thermal-equilibrium LGM and CMD 



models. To this end, the temperature dependence of 
the mean yield (< n >)of individual clusters is investi- 
gated. In order to compare a recent well-known thermal 
Arrhenius-type plot in nuclear multifragmentation phe- 
nomenon H ||, we plot In < n > versus 1/T. Fig. || gives 
a family of these plots for the disassembly of 129 Xe within 
the framework of I-LGM (left panels) and I-CMD with 
Coulomb interaction (right panels). Again, as Fig. [l], the 
clusters are classified according to their masses (upper 
panels), the light isotopes (the middle panels) and the 
charge numbers (the lower panels). For all the panels, 
the obtained Arrhcnius plots are linear for the lower T 
side, and their slopes generally increase with increasing 
A or Z value. Generally, the thermal scaling is expected 
when the yields, for a fixed nucleon number system, are 
dominated by fragment binding. This is the case when 
the temperatures are low compared to the binding energy 
per particle. At these temperatures, one can anticipate 
one large fragment surrounded by many small clusters. 
However, contrary tendency reveals in the high T side 
where In < n > increases with ie. decreases with 

increasing T. In this case, nuclear Arrhenius plots of 

< n > with 1/T are not valid but the Poissonian re- 
ducibility still remains (see Fig. |^). This behavior of 

< n > at higher T is related to the branch of the fall 
of the multiplicity of IMF (Nimf) with T where the 
disassembling system is in vaporization |2^, M and 
hence only the lightest clusters are dominated and the 
heavier clusters become fewer and fewer with increasing 
T. The temperature where the Arrenhius-type plot be- 
gins to deviate from the linearity just indicates the onset 
of transition from the Fermi liquid phase to gas phase in 
I-LGM and I-CMD Q. In other words, the Arrenhius 
law looks valid only below the critical phase transition 
point, ie. in the nuclear liquid phase. Recognizing this 
phenomenon, in the following sections we only focus on 
the branch of lower temperature (ie. the liquid phase) 
where the thermal scaling exists to discuss the Arrhenius 
law and their slopes. 



C. "Emission Barriers" 

1. Model dependence 

From Fig. |3] the slope parameter can be directly ex- 
tracted in the lower T side as a function of Z or A. In Ref. 
H Moretto et al. has interpreted these slope parameters 
as "emission barriers" of specific individual fragments. 
Fig. |] gives the emission barrier of individual fragments 
with different A, ISO and Z in the framework of I-LGM, 
I-CMD with/without Coulomb interaction. The error 
bar in the figure represents the error in the extraction of 
the slope parameter. The first indication from this fig- 
ure is that the emission barrier in the I-LGM case is the 
nearly same as the I-CMD case without Coulomb force, 
which supports that I-LGM is equivalent to I-CMD with- 
out Coulomb interaction rather well when the nuclear 
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FIG. 3: Arrcnhius-type plot: the average yield per event of differ- 
ent clusters classified with A (top), ISO (middle) and Z (bottom) 
as a function of 1/T. The left panels show the results with I-LGM 
calculation and the right present the results with I-CMD including 
Coulomb force. The solid lines are fits to the calculations using a 
Boltzmann factor for < >. The symbols are illustrated on the 
figure. 



FIG. 4: The emission barrier extracted from the Arrhenius-typc 
plots as a function of cluster mass (top), isotopic mass (middle) 
or cluster charge (bottom) in the cases of I-LGM (solid squares), 
I-CMD without Coulomb (solid circles) and with Coulomb (open 
circles). The unit of the emission barrier is MeV throughout this 
paper. The solid lines are fits with the Eq. (8) or (9), and the dot- 
dashed lines represent the fits with the Eq. (10) or (11). 



potential parameter is moderately chosen, but I-LGM is 
a quick model to analyze the behavior of nuclear dissoci- 
ations. The inclusion of long-range Coulomb interaction 
makes the emission barrier of individual fragments much 
lower since the repulsion of Coulomb force reduces the 
attractive role of potential and hence make clusters es- 
cape easily. The second indication is that the emission 
barriers increase with A (Z) at low A (Z) values and tend 
to be saturated at high A (Z) ones. Similar experimental 
results have been observed for individual fragments with 
different Z in Ref.||] or different A in Ref.[Q. However, 
the middle panel of Fig. 4 shows that bare dependence of 
emission barrier of ISO on A in the fixed atomic number 
Z, which indicates that the Z dependence of barrier is 
perhaps more intrinsic and the A dependence is basically 
due to the average effect over the species with the same 
A but different Z. 



2. Source size and Coulomb interaction dependence 

On the origin of these barrier, the surface energy and 
Coulomb energy would play the roles. If the cluster emis- 
sion is mainly controlled by its surface energy, it would 
suggest barriers proportional to Z 2 / 3 (A 2 / 3 ). In the case 
of I-LGM and I-CMD without Coulomb, we can try to fit 
the barrier for the particles with different mass number 
by 

Bcoui.off ^^x Ai 2/3 , (8) 
or for the particles with different charge number by 

Bcoui.off = ci x ((A/Z) ftt * Z,) 2 ' 3 , (9) 

where (A/Z) fu is a fit coefficient of A/Z for emitted par- 
ticles, and Ai (Zi) is the mass (charge) of particle. c\ is 
the fit constant for surface energy term. The solid line 
in the Fig. ^a is a function of Eq.(8) with c\ = 8.469 and 
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the solid line in the Fig. [i|c is a function of Eq.(9) with C\ 
= 8.469 and (A/Z) fit = 1.866. These excellent fits imply 
that the surface energy play a major role in controlling 
the cluster emission when the long range Coulomb force 

I 



for the particles classified with different mass number, or 
I 



for the particles classified with different charge number, 
where c 2 is a fit constant for surface term and rc ou l is 
chosen to be 1.22 fm. Z res is a fitted average charge 
number of the residue. (A/Z)fn is chosen to be 1.866, as 
taken from the fits for I-LGM. The overall fits for A and 
Z dependent barrier in the case of I-CMD with Coulomb 
force give c 2 = 12.921 and Z res ~ 41 with the dot-dashed 
line in Fig. 4a and 4c. The excellent fit supports that 
the Coulomb energy plays another important role in the 
cluster emission. 



In the case of I-LGM and I-CMD without Coulomb, 
one would expect the barrier for each Z (A) to be nearly 
independent of the system studied if only the surface en- 
ergy is substantial to the emission barrier. The left panel 
of the Fig. U shows the results for Ba, Bjso and Bz for 
three different systems in the I-LGM case. The same 
freeze-out density of 0.38po and the same N/Z is chosen 
for the systems of A sys = 80 and A sys — 274. Actually, 
it appears to have no obvious dependence of emission 
barrier on source size as expected for the role of surface 
energy. The solid line in the figure is the same as in 
Fig. 0. However, when the long-range Coulomb interac- 
tion is considered, the emission barrier reveals a source 
size dependence. The right panel of Fig. ^| gives the emis- 
sion barriers Ba, B/so and Bz in the case of I-CMD with 
Coulomb force. It looks that the barrier increase with the 
decreasing of charge of system, which can be explained 
with the Eq. (10) and (11) where the decreasing of the 
residue Z res will result in the decreasing of the Coulomb 
barrier and hence the increasing of the emission barrier. 
The lines represent the fits with the Eq.(10) and (11) for 
three different mass systems. 



is not considered. However for the cluster emission with 
the Coulomb field, we can assumed that the barrier is 
mainly constituted by the surface energy term and an 
additional Coulomb term as 



(10) 

I 



(11) 

I 

3. "Freeze-out" density dependence 

In the above calculations, the freeze-out density of sys- 
tems is fixed at ~ 0.38po- Considering the freeze-out den- 
sity is an important debating variable in the latter stage 
of heavy ion collisions, here we will discuss the possible 
influence of freeze-out density on the emission barrier of 
clusters. The calculations at the freeze-out density of 
0.177p and 0.597p for 129 Ae, corresponding to 9 3 and 
6 3 cubic lattices respectively, are supplemented to com- 
pare. Fig. U gives the results of Ba, Bjso and Bz at dif- 
ferent density. It looks that there are no obvious freeze- 
out density dependence in the both cases of I-LGM and 
I-CMD. This is also consistent with that assumption that 
the surface energy is the dominant role in controlling the 
cluster emission. 



IV. CONCLUSION 

In conclusion, the Poisson rcducibility and thermal 
scaling of the emitted clusters is explored in the lat- 
tice gas model and molecular dynamical model. It in- 
dicates of a statistical nature of such cluster emission. 
The Poisson reducibility illustrates that the clusters are 
produced independently each other and stochastic. But 
near the liquid gas phase transition, the large fluctua- 
tion breaks down the Poisson reducibility. The ther- 
mal scaling is existed when the temperatures are low 
compared to the binding energy per particle. At these 
temperatures, one can anticipate one large fragment sur- 
rounded by many small clusters, ie. the nuclei is in the 
liquid phase. The calculations looks qualitatively consis- 
tent with the recent experimental observation of Poisson 
reducibility and thermal scaling by Moretto/Wozniak's 



R v 4 2/3 1.44 x Ail {A/ Z) fit x Z re 

Bcoul.on — C2 X A; 



rcouiiAi 1 / 3 + ((A/Z) fit * Z res y/*) 



Bcoul.cn = C 2 X {{A/Z) fu * Z t f IZ 



1.44 x Zi x Z r 



r C oui((Zi * (A/Z) flt y/z + (Z res * (A/Z) flt y/3) 



7 




FIG. 5: The source size dependence of the emission barriers for the 
different clusters classified with mass (top), isotopic mass (middle) 
or cluster charge (bottom) from in the cases of I-LGM (left panel), 
I-CMD with Coulomb (right panel). The lines in the left panel are 
fits with the Eq. (8) or (9), and the solid, dot-dashed and dotted 
line in the right panel represents the fits to the emission barrier of 
A sys = 80, 129 and 274, respectively, with the Eq. (10) or (11). 



FIG. 6: The emission barrier of 129 Xe for the different clusters 
classified with mass (top), isotopic mass (middle) or cluster charge 
(bottom) at the different freeze-out density in the cases of I-LGM 
(left panel), I-CMD with Coulomb (right panel). The lines arc fits 
with the Eq. (8) or (9) in the left panel and with the Eq. (10) or 
(11 ) in the right panel. 



group even though the system studied is different and 
the temperature was supposed to be proportional to the 
total transverse energy in their experiments. This also 
supports somehow that the lattice gas model and classi- 
cal molecular dynamics is a useful tool to simulate the 
nuclear disassembly. 

Further, based on the Arrhenius law in the liquid 
phase, we extracted the emission barrier for clusters with 
the different mass, light isotope mass and charge num- 
ber. Also, the systematic fits with the formula embodied 
with the surface energy and Coulomb interaction were 
performed and the overall good fits were reached. The 
results suggest that the cluster emission is mainly con- 
trolled by both the surface energy and the Coulomb in- 
teraction. In the framework of the lattice gas model and 
molecular dynamics model without the Coulomb inter- 
action, the emission barrier relies on the cluster charge 
with the Z 2 / z {A 2 / z ) law and it does not depend on the 
source size and freeze-out density, which indicates that 
the surface energy play a basic dominant role to con- 



molecular dynamics model with the Coulomb force, the 
emission barrier will decrease strongly according to the 
Eq.(10) and (11) and it decreases with the increasing of 
the source size, illustrating that the Coulomb interac- 
tion also play another weighty role to control the cluster 
emission. 
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